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THEORETICAL LIFT AND DRAG OF THIN TRIANGULAR WINGS AT SUPERSONIC SPEEDS 

By Clinton E. Brown 
December 19 b 6 


Equation (29) of the section "Drag Due to Lift" (pp. 10-12) has been 
found to be in error as a result of the omission of a compressibility 
factor in the treatment of the two-dimensional leading- edge suction force. 
This factor can be obtained by application of the well-known Prandtl— 
Glauert rule by which incompressible solutions can be corrected to obtain 
the corresponding compressible— flow solution. For the two-dimensional 
case it is found that the effect of compressibility on a flow having a 
given vcrticity distribution in a plane is to reduce the velocities nor- 


mal to the surface by the factor ,J 1 — M* 2 where M* is the Mach num- 
ber of the flew. Therefore, if the strength of each vortex element is 


increased by a factor . ~z , the resultant velocities normal to 

•J 1 - M* 2 

the surface will again be equal to those of the incompressible flow. In 
this case, however, the tangential velocities at the surface and therefore 


the forces on the surface are increased by the ratio — rr:: . This 

■J 1 - M* 2 

concept is well known in thin— airfoil theory where it gives the result 


that the lift— curve slopes of thin airfoils are increased hy 

STi- 

and also that the total resistance for M < 1 is zero. This result also 
indicates that the leading— edge suction force has "been increased over 

that of the incompressible flow by the ratio — # 

•J 1 - M‘ 2 


It appears that the incompressible equations governing the leading- 
edge suction force (equations (21) and (22)) must be corrected as follows: 


The leading-edge suction force on a two-dimensional plate will be 


F n 


PnG 2 
J 1 - M* 2 


( 22 ) 


when the vorticity distribution 7 at the leading edge or the tangential 
velocities are given by the following relations: 


n 


7 

-> 0 


2G 


J 1 - M’V s n 


- 2 - 


or 

G_ 

8 n /. 1~ M ,S */e n (21) 

If the corrected relations are now carried through the analysis as given 
on pages 10, 11, and 12, the final result for the drag coefficient Cj^ 

given in equation (29) should read 

- vs K 1 *!)- _ 

In this case M* is the Mach number component normal to the leading 
edge and can he shown to he 


“5 J l - 3 2 c 2 

J 1 + C2 

Equation (29) therefore becomes 



It should he noted that, because of the infinite pressure leading 
at the tip of the triangular wing, the suction force indicated will 
probably be nonexistent in the very highly leaded tip regions, as is 
the case at subsonic speeds with sharp— nose airfoils at finite angles 
of attack. The amount of suction force obtainable in practice will de- 
pend on the leading-edge radius and angle of attack and, in the final 
case, must be obtained experimentally. 
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THEORETICAL LIFT AND DRAG OF THIN TRIANGULAR 
WINGS AT SUPERSONIC SPEEDS 
By Cl in ton E . Brown 


SUMMARY 


A method is derived for calculating the lift and the drag due 
to lift of point-forward triangular wings and a restricted series 
of swept-back wings at supersonic speeds . The elementary or 
supersonic source solution of the linearized equation of motion 
is used to find the potential function of a line of doublets. The 
flow about the triangular flat plate is then obtained by a surface 
distribution of these doublet lines. The lift-curve slope of 
triangular wings is found to be a function of the ratio of the 
tangent of the apex angle to the tangent of the Mach angle . As the 
apex angle approaches and becomes greater than the Mach angle, the 
lift coefficient of the triangular wing becomes equal to that of a 
two-dimensional supersonic airfoil at the same Mach number. 

The drag coefficient due to lift of triangular wings with 
leading edges behind the Mach cone is shown to be close to that of 
elliptically loaded wings of the same aspect ratio in subsonic 
flight. The resultant force on wings with leading edges outside the 
Mach cone, however, is shown to act normal to the surface and thus 
an induced drag equal to the lift times the angle of attack is obtained. 


INTRODUCTION 


In reference 1, Jones calculated the lift of thin point-forward 
triangular wings for the cases in which the apex angle of the wing 
was very small. It was pointed out that the results obtained should 
be applicable in both supersonic and. subsonic flight, the criterion 
for the case of supersonic flight being that the apex angle be small 
as compered with the Mach angle of the flow. The present paper, 
making use of less restricted theory, extends Jones' work to the case 
of triangular wings having large apex angles and traveling at super- 
sonic speeds. A recent paper was published by H. J. Stewart 
(reference 2) in which the lift of triangular wings has been computed 
but the method used appears to be entirely different. 
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In the present theory, the linearized equation of motion was 
used and the results must therefore he restricted to small angles 
of attack and moderate supersonic Mach numbers . The solution 
which has been found should hold good for large values of the apex 
angle up to and coincident with the Mach angle. Jones (reference 3) 
and Puckett (reference 4) have found solutions for the drag of 
triangular wings of small thickness at zero angle of attack. The 
solutions are applicable to wings having the leading edges either 
in or out of the Mach cone springing from the apex of the wing. 
Puckett has pointed out that, for the case where the leading edge 
is ahead of the Mach cone, these solutions can also be used to 
calculate the lift; thus, with the present solution, the lift for 
the whole range of ap9x angles at supersonic speeds may be obtained. 
The pressure distributions and lift-curve slopes obtained in the 
present paper can be used to obtain the lift and drag characteristics 
of a limited series of swept-back wings. The drag due to lift of 
the triangular wing has been calculated and a suction force has 
been found to exist on the leading edge . In order to use the suction 
force, however, it appears necessary tc provide an airfoil section 
with a rounded leading edge . The aut.. r is indebted to 
Mr. Arthur Kantrowitz of the Langley Memorial Aeronautical Laboratory 
for suggesting the method used to calculate the induced drag. 


a 


angle of attack 



b maximum span of wing 


p '/M 2 - 1 


C tangent of apex angle 


C. 


'L 


lift coefficient 



Cp drag coefficient due to lift 



Dj drag force due to lift 
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E source strength 

c apex angle of wing measured, from flight direction 
f(o) double t-distribution function 
F suction force on wing leading edge 



i = - 1 


strength of line dov.blet 


a = — 


U 


length of wing or root chord 
lift force 


X = 

M 

M- 

Ap 

q 

E 

p 

s 

s n 

u 


(./_ 




S 2 c 2 


tanh 


c '/c 2 - a 2 '/l - 3 2 a 2 

Mach number 
Mach angle ( sin 


\/l - 3 2 cj 2 da 


( sln_1 s) 


lifting pressure 
dynamic pressure 


it' 1 *) 

distance along leading edge from wing at>ex 
density of the fluid 
wing area 

distance normal to leading edge 
velocity increment in x-direction 

WJ 

velocity increment normal to leading edge / 

V^n/ 


k 
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v 


V 

w 


W 

x, y, z 
X 1 z x 


0 

<*o 

02 


velocity increment in y-direction 



flight velocity 

velocity in the z-direction (£-&) 

resultant velocity in z-direction created "by 
the doublet distribution 

coordinates of an arbitrary field point 

coordinates of a source or doublet 


disturbance-potential function 
potential of a supersonic source 
potential of a line of supersonic sources 
potential of a line of supersonic doublets 


Subscripts : 

n normal to leading edge 

A triangular-wing condition 

od infinite-span or two-dimensional wing condition 


THEORY FOR LIFTING TRIANGLE 

The linearized equation of motion of a nonviscous compressible 
fluid may be written 


+ <£$ _ (m2 - 1 ) <£& 
ciy 2 dz 2 Sr 2 


where 0 is the potential function assumed to represent the effect 
of a small disturbance set up by the body being considered. The 
body in this case is a triangular flat plate having its vertex at 
the center of the coordinate system and lying in the xy-plane (fig. 1) . 
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The problem is to find a solution of equation (l) that will satisfy 
the known toundary conditions which are (l) that the flow he 
quiescent ahead of the Mach cone and (2) that the flow at the surface 
of the plate he tangent to that surface . Because of the linear 
character of the differential equation (l), more general potentials 
can he but It up from simple well-known solutions such as that for 
a single source 


0 O . — ■ (2) 

+ T) 


where p = |/ M 2 - 1# The potential of a line of sources with 
strength proportional to x can he found as follows: 


r Ex, dx, 

01 = j ... (3) 

u 0 \/fx - x-l)- - p rJ (y - ax x ) ' - |3 2 (z - Tx-jJ 2 
y z 

where a = — ^ t = — , and x' is the value of x n for which 

X 1 X 1 

the denominator of the integrand is zero. Physically interpreted, 
the range of integration is from the origin to the last source point 
which can influence the field point. Performing the integration 
yields 



v- 


( 


p p 

x - P gy - a T Z 

1 - - PV) 372 


ctnh 


-1 


x - |3 2 q.y - P 2 tz 


(1 - |3 2 a 2 - |3 2 t 2 ) x 2 - P 2 (y 2 + z 2 ) 


GO 


If two such source lines of opposite strength ere brought together 
from the z-direction at the xy-plane while the product of source 
strength and the angle between them is kept constant, the potential 


V" 


6 
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of a line of doublets in the xy -plane at an angle tan'^a from the 
x-axis is obtained. Thus differentiating with respect to T and 
setting T = 0 gives 


$2 



ctnh"' 1 * £ 


(5) 


where 

t = - P 2 °.y 

/ (1 - 3 2 <t 2 )|x 2 - 3 2 (y 2 + z 2 )J 


and I is the doublet strength. Differentiating the potential 
function with respect to z gives the vertical velocity w: 


v 




_ i£ 

- P 2 a 2)3/2 



ctnh - -^ 


+ 2iz 2 — P_ 2 .L-<; ..- 3 2 o.y) \/y 2 - P 2 (y_ 2 + z 2 ) 

'| (x “ 3"ay) 2 - jjL - 3 2 cr2j j ^2 _ p 2 (y2 + z 2)| |> 
s ’* \J 


( 6 ) 


It will be noticed that the line doublet creates a conical field as 
the velocity is only a function of z/x and y/x. Since the 
triangular -flat plate is a conical body which creates a conical field, 
an attempt will be made to build up the flow about the lifting 
triangle by a suitable distribution of line doublets inasmuch as the 
addition of two or more conical fields having the same vertex always 
creates another conical field. The distribution of line doublets 
must satisfy the boundary conditions at the body surface which may 
be written: 


W = Va 
z=0 


( 7 ) 
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where W is the resultant vertical velocity of the line-doublet 
distribution. If the distribution of line doublets is f(cr) 

V = f (a) w (a; 1; z\ da (-C < - < c) (3) 

z— ? 0 J- c V X 3C/ V x ' 

in which tan "C = S, the angle of the leading edge. The distribution 
function f (a) can be found in a rather simple way by analogy with 
the solution for incompressible two-dimensional flow about a flat 
plate. Differentiating equation (8) with respect to y/x and 
setting z = 0 gives 



which i3 of the same form as the integral equation obtained when the 
incompressible flow normal to a two-dimensional flat plate is 
constructed by a doublet distribution . The expression f(a) for 
the incompressible case would be 



( 10 ) 


That this expression is a solution of equation (8) must be verified 
by substitution in equation (8), inasmuch as equation (9) is a 
divergent integral. This proof is carried out in appendix A. The 
value of the velocity in the x-direction u can now be obtained 


f(a) °^1 da (11) 

3x 

The integration indicated in equation (ll) is presented in appendix B. 
The expression obtained from equation (ll) on the lifting surface 
(z = 0) gives 




u 


( 12 ) 
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the sign of the expression being opposite for the two sides of the 
plate. The result in equation (12) shows that the shape of the 
pressure distribution is independent of aspect ratio. 

The lift for the isosceles triangle with root chord c r is 


= / Ap dS = P 
1 / . / 


pVuCy. dc 


V-C 


( 13 ) 


and substituting equation (12) in equation ( 13 ) and integrating gives 


L = pV Ic r 2 jt 2 C 2 


(14) 


and 


, L 2ljt & C 

X " qS V 


( 15 ) 


The value of the constant I must be obtained by solving equation (8) . 

The value of the normal velocity at the plate and hence the 
angle of attack may be found by integration of equation (8) and 
letting z = 0. The integration is involved and the method of 
integration is given in appendix A. The resulting expression for W 
is obtained as 



R 2 2 
. JtuSL - 


a® /] 


tanh' 


- 1 \/ 


pV 



( 16 ) 


From this equation the value of I may be calculated. If the numerator 
and denominator of the integrand are multiplied by p the resulting 
integral can be seen to be dependent upon only the quantity pC 
or tan tj tan |i. The value of the integral may be obtained easily by 
making the substitution P^c2 - p 2o2 = n 2 ani plotting the resultant 
expression. This procedure has been followed for values of PC 
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between 0 and 1 and the result Is given in figure 2. The 
value of I is found to be 

I = (17) 

Tt + X 

where X is the integral term of equation (l 6 ) . The lift-curve 
slope is now 


£l _ 2it p C 
a Jt + X 


(18) 


As (3C approaches zero, X also approaches zero and the lift-curve 
slope from equation ( 18 ) approaches' Jones 1 value. (Reference 1.) 
Equation (l 8 ) shows that the lift-curve slope is a function of only 
the apex angle and the parameter tan e/tan g . It is interesting 
to note that mathematically there is a finite lift-curve slope at 
the Mach number 1.0. The theory is not valid, however, near M = 1 
because of the original assumptions used in obtaining equation (l) . 
The lift-curve slopes of two triangular wings are plotted in 
figure 3 against Mach number. 

Rewriting equation (l 8 ) in a convenient form gives 


/ c T \ 2 n 2 

J _L\ = tan u 

v a /A (it + x) \/m 2 - 1 


(19) 


Now, according to Ackeret's result (reference 5) the lift-curve 
slope of a two-dimensional flat -elate is 



A single curve for all Mach numbers can therefore be plotted if the 
ratio tan c/tan g is used for the abscissa and ^/aj 

is used as the ordinate . This curve is shown in figure 4 . 


10 
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It can "be seen that as the anex angle approaches the Mach angle 
the triangular wing provides the same lift coefficient as a two- 
dimensional wing at the same Mach number. 

The case of the triangular wing having the leading edge ahead 
of the Mach cone from the apex has been treated in reference k. 

It was found that the lift coefficient obtained is the same as that 
of a two-dimensional airfoil flying at the same Mach number. The 
curve shown in figure 4 therefore becomes flat at values of 

> 1 . A typical pressure distribution over a wing having 

fo— - >1 is shown in figure 5. 
tan u 


The thin-airfoil theory used herein gives the result that the 
resultant force is directed normal to the plate, a result quite 
like that obtained from the thin-airfoil theory at subsonic speeds. 
In the solution for subsonic speeds, however, a simple extension 
described in reference 6 permits calculating the force due to the 
suction on the leading edge. It is reasonable to suppose that the 
same method is feasible for the triangular wing, as the pressure 
distribution in the neighborhood of the leading edge is identical 
in the limjt with that for a two-dimensional flat plate in subsonic 
flow; that is, the velocities nornal to the leading edge are in the 
form 


-where s n is the distance measured normal to the edge and G is a 
constant. According to reference 6 the force normal to the edge in 
the direction of the velocity is 


The value of the velocity in the y-direction on the triangular wing 
has been calculated to be 


BRAG BUE TO LET 


s n ~» 0 \A n 


( 21 ) 



( 22 ) 


In Z 


X 



V 
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Combining this expression with equation (.12) gives 



(?3) 


where R - x Vl + C 1 " . The normal force on a small element dR of 
leading edge "becomes 


c 9 

n _ P rtCBr rr '- 
<3R 2 




and the force on one edge of the isosceles triangle with root 
chord c r is : 

r \A + c 2 

EcLR 



P it 3 I 2 Cc r 2 (l + C 2 ) 
- 


(25) 


and the force in the flight direction from two edges "becomes 


F = 2F n sin e 


pir 


3n2 


V 1 


+ C 2 c^ 2 ! 2 


2 


( 26 ) 
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Substituting from equation (l4) 



La A / 1 + C 2 \ 



( 27 ) 


vhere b is the maximum span of the triangular wing. The induced 
drag or, more exactly, the drag resulting from the lift may be 
written 


It will be noticed that this result is identical with Jones 1 result 
(reference 1) in the limiting case of C = 0. The induced drag 
coefficient is found to be 


where A is the aspect ratio. Equation (29) indicates that the 
triangular wings can obtain a considerable suction force at the 
leading edge and that the drag coefficient due to lift of th 3 wings 
is very close to that obtained from elliptical wings of the same 
aspect ratio at subsonic speeds. It should be pointed out, however, 
that as soon as the wing leading edge passes through the Mach cone, 
the possibility of obtaining a leading _ edge suction is gone and the 
resultant force must become normal to the plate surface. This 
transition corresponds quite similarly to the case of a two-dimensional 
airfoil passing through the speed of sound. 


D i = La - F 
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DISCUSSION AND CONCLUSIONS 


The lift at supersonic speeds of triangular wings having 
straight trailing edges has teen shown to approach the lift of a 
two-dimensional airfoil as the leading edge approaches the Mach cone 
scringing from the apex of the triangle . For the case where the 
trianguler wing lies "behind the Mach cone, a suction has been found 
to exist on the leading edge . In order to utilize this suction 
force in practice it would appear necessary, as in subsonic flow, 
to provide an airfoil section vi th a rounded leading edge. Triangular 
wings ehould be capable of higher L/D ratios than unswept wings 
at supersonic speeds when operating with their leading edges not too 
far behind the Mach cone; the improvement should be due to both 
reduced wave drag and reduced induced drag. 

The lift and drag of a series of limited swept-back wings 
may also be calculated with the method developed. It will be noted 
that the pressure distribution over the triangular wing cannot be 
changed if the trailing edge is cut off from the tip to the center 
line along an angle always greater than the Mach angle . This fact 
arises from the nature of the supersonic flow in which disturbances 
cannot propagate any farther forward than the Mach cone from the 
origin of the disturbance. The aforementioned procedure produces 
therefore a series of tapered swept-back wings having pointed tips . 

A new series can also be constructed by cutting off the tips along 
lines having angles greater than the Mach angle. In each case the 
pressures over the remaining portions of the wing will be the same 
as though the cutbacks had not been made . 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., Nov. 29, 1946 
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APPENDIX A 


CALCULATION OP VERTICAL VELOCITY INCREMENT 


The value of W is, from equations (6) and (8), 

p* p 

w » I / vfc2 - a 2 / 1 £-22— j-lotah-lf - _J 

J, Li k 2 -d 

V t 


+ 2z 2 do - .P 2 .(.y 2 + ^ 2 ) 

[a 2 -'* 2 - p 2 z 2 ) - 2xya + y 2 + z^j 2 


(Al) 


Integrating by parts gives 


W « I 


< 


■j£s-- L 
,e ' 


./c 2 - 


t> (1 - p 2 


a 2 ) 3/2 |I 


P- 

Ls . 


- ctnh”\f 


da > do 


+ I \j£- p 2 (y 2 + z 2 ) 


f- 


-flSJlJL 


-C 


\/c 2 - o 2 I 


d2 2\ , 2 , 2 

p z ; - cxya + y + z 


L 


x 


tan 


n n n 

-1 ( X ~ p T ).T ~ 


do 


(A2) 


z’y x 11 - p 2 (y 2 + z 2 ) 


x 2 - f>V * 


Z 2 \ 
L i 


The .integration of the term under the indefinite integral can he performed by parts to 
obtain the result 


. — 2 (— L_ 

'/: - ru s \i 2 - 1 


- ctnh -1 £ ' 


\ih£ - B g (y 2 + ■> 'r:zh- 2 ' - 26 2 / 2 •• f>, 2 s 2 y + f^.Wy 2 t JT“_)_ 

|3 - •, x 2 - p 2 z 2 ; | o' - iff 1 ) - 2' ye + y" + z" 

I 


... p p. 

+ +on“l 3 ^ -‘) cr - X 7 


P 2 z 


tan 


z V x 2 - (3 2 (y2 + ) 

Substituting in equation (A2) and rearranging terms gives: 
/oC 


V = I 


o 2 

a- 6 


.. /ctnh“^-f - -- - ~ 

■C >/l - (>V '/c 2 - c 2 y 


t 2 -l 


da 


+ iiA 2 - p^-i - + z*) 

U 2 - |b 2 Z 2 ) 


' J2- ^ „2> 


/V 


-C 


Vo 2 - a 2 


2 ) 2[x 2 - p C (,v 2 + z^llax - yfx 2 - ft 2 y 2 - 2g 2 z 2 ) 

a 2 (x 2 - p 2 z^) - 2xya + y 2 + z^ 


dcr 


(A3) 


VJ1 
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It vill be noticed in the preceding operation that all the terms 
containing the singularity of the form l/z cancel. If z is 
made to approach zero and terms are collected, equation (A'i) 
becomes 


I 



Completing the integration of the second terms gives 



Differentiating expression (A5) with respect to (y/x) and performing 
the integration gives 



the vertical velocity is therefore constant over the plate surface and 

the expression f(a) = v/c^ - 0 ^ is truly a solution of equation (8). 
It is possible that this solution is not unique; however, other 
solutions would undoubtedly lead to physically impossible conditions . 
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APPENDIX B 

CALCULATION OF THE AXIAL VELOCITY INCREMENT 


The value of the x-component of velocity may "be written from 
equation (ll) 


u 




2 2 

ovx - y~ - z 


(a 2 (x 2 -p^ 2 ) - 2cyx + y 2 + zf} 


■=71 d0 


(Bl) 


This integral can he broken up into two separate integrals, as follows: 


u = 


Vx 2 - b + z 2 ) 

(x 2 - 


xy 


r "' C \/c 2 - c 2 (?-g - 


_ - 

x 2 - j 5 ?, 


-C 


-^L g — + _£ 


2 , 2 \2 

+ z \ 


da 


\ 


x 2 - p 2 z 2 x 2 - P 2 Z 2 . 


y 2 + z 2 


p p 


r 


X 2 - p^z 2 ; 


-C ^ 


Vc 2 - a 2 da 


Txycr 


X 2 - P 2 Z 2 


2 2 \ 2 
y + z 

x 2 - P 2 - 2 


(B2) 


The first integral appearing in equation (B2) can he integrated hy 
parts to give, for the first complete term, the integral 


- Izxy y^ 2 - P 2 (y 2 + z 2 ) 

(j? - e 2 , 2 ) 2 




dtr — ^ 

2xyq + y 2 + z 2 \ 

. p 2 z 2 ^ . p 2 z 2 J 


(B3) 
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Evaluation of this integral, which may he found in reference 7 , 
equation (228), gives finally 


j. IXVrt J 

1 xy - 

iz,/ x 2 - g 2 (.y 2 + z 2 ) 


2(3? - | 3 2 z 2 ) N 

p- 

xv - i z V' x 2 - B 2 (y 2 + z 2 ) 

1 

2 


( 

0 co 

x^ - p 2 z 2 



V. 




+ - P 2 (y 2 + z 2 ) 


1 / 


'c 2 - 


xv 


+ iz v /x 2 ~ 3 2 (v‘~ + zj 1 

x 2 - P 2 Z 2 


> 


(B4) 


the sign of the expression being opposite for the two sides of the 
"olate . The second integral term of equation (B 2 ) can "be integrated 
by breaking the integrand up into four partial fractions, as follows: 


vte - a 2 


(x 2 - ( 3 2 z 2 )" da 


\kiz3\/lx 2 - p 2 (y 2 + z^)j J 


a - 


l 


XV - iz’\/x 2 - p 2 (y 2 + z 2 ) 


x 2 - ^z c 


(x 2 - g 2 ? 2 ) 1 " da 


4 £.J\£ - 0 2 (y 2 + z e )f 


0 - xv - iz l/:^ - $ {y d + z 2 ) 
x 2 - P 2 z 2 


(x 2 - g^z 2 ) 2 da 


4 iz 3^1 j? - p 2 (y 2 + .z 2 )] 3 


xv + iz/x 2 B 2 (.y 2 +_.z^). 


2 .0 2 

XT - ( 3 '-z 


(x 2 - B 2 z 2 ') 2 da 


4z 2 | / /jx 2 - P 2 (y 2 + z 2 )J 


a - xv + Izi/x 2 - B'-fv 2 + z L ) 
x 2 - 3 2 z 2 


05 ) 


2 \ 
/ 
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The expressions (B5) may now he integrated (reference 7> equation 207) 
giving the expression for the complete second term of equation (B2) 


+ si < U - 


2 '/ 


xy - iz '/x 2 ~ ,6^ (y 2 + z'") 


2 


x 2 - P^z 2 


/ 




+ iz V'x 2 ~ 6 2 f y2 + z 2 ^ 


X 2 - pV 


+ i'-'Tz V' 7 ? 2 ~ P 2 (y 2 + z*" - ) ! xy - iz'/x 2 - P 2 (.y^ + z^) 
2(x 2 - p 2 z 2 ) 2 ' / 


;A 2 


xv - Iz y/x 2 - P^(y 2 + z f ~) 


- p 2 z 2 


xy + iz\/ x 2 - B^(y 2 + z 2 ) 


(B6) 


i2 . 


V 


xy + iz vr - (y 2 + z 2 ) 

‘ x 2 - p 2 z 2 


Combining equations (B^) and (b 6) and setting z = 0 yields for u 
on the surface 



u 


(B7) 


) 
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Figure 1.— Coordinate system. 
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Figure 2.- Variation of A with H2JL 
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Figure 3.- Variation with Mach number of Cl/cc for two triangular wings. 
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Figure 5.- Typical pressure distribution for triangular wings with 
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